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ABSTRACT 

We investigate, independently of specific emission models, the constraints on the value 
of the bulk Lorentz factor V of a fireball. We assume that the burst emission comes from 
internal shocks in a region transparent to Thomson scattering and before deceleration 
due to the swept up external matter is effective. We consider the role of Compton drag 
in decelerating fast moving shells before they interact with slower ones, thus limiting 
the possible differences in bulk Lorentz factor of shells. Tighter constraints on the 
possible range of T are derived by requiring that the internal shocks transform more 
than a few per cent of the bulk energy into radiation. Efficient bursts may require a 
hierarchical scenario, where a shell undergoes multiple interactions with other shells. 
We conclude that fireballs with average Lorentz factors larger than 1000 are unlikely 
to give rise to the observed bursts. 
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1 INTRODUCTION 



In the last few years, increasing evidence in favor of the fire- 
ball model (Rees & Meszaros 1992) for gamma-ray bursts 
has been gathered, thanks to the observations of BeppoSAX 
(Boella et al. 1997). The power law decay of the optical 
afterglow of several bursts has lasted for time-scales of a 
year, in striking agreement with the simplest fireball sce- 
nario (Wijers, Rees & Meszaros 1997), in which a shock 
wave propagates in the interstellar medium (ISM), acceler- 
ating particles which then emit by the synchrotron process 
(Sari, Piran & Narayan 1998). However a single shock wave 
cannot account for both the temporal behavior of the 7-ray 
emission (Fenimore et al. 1999a) and for the requirement of a 
high efficiency in converting its kinetic energy into radiation 
(Sari & Piran 1997). In fact the time variability structure of 
bursts appears to be constant during the whole of the high 
energy emission (Fenimore et al. 1999b), while the deceler- 
ation of a shock slowed down by interactions with the ISM 
would produce a time dilation between the first and the last 
spike of the burst. Moreover the observed variability cannot 
be explained as the consequence of inhomogencities in the 
ISM since a radiative efficiency of less than 1% would be 
expected (Sari & Piran 1997; see however Dermer, Bottcher 
& Chiang 1999). 

In the internal shock scenario, put forward by Rees & 
Meszaros (1994), the inner engine produces many relativis- 



tic expanding shells (or an unsteady wind) with a distribu- 
tion of Lorentz factors centered on a mean value (F) and 
width AF ~ (r). The burst radiation is produced through 
the dissipation occurring when a faster shell catches up a 
slower one. No matter the physical mechanism producing 
the observed photons, the inner (hidden) engine is responsi- 
ble for the temporal structure of the observed burst through 
the time history of the shell emission (Kobayashi, Piran & 
Sari 1997). Even if we may have information on the time- 
scale of the shell ejection, it is hard to estimate the bulk 
Lorentz factor (F) of the relativistic outflow. This is how- 
ever a critical parameter to unveil the mechanism(s) that 
powers the outflow and the radiating process that produces 
7-ray photons. In fact, the knowledge of (F) allows us to es- 
timate the amount of baryon loading of the fireball and the 
intrinsic frequency of the emitted photons. The constraints 
related to the compactness problem (see e.g. Piran 1992) 
give a lower limit (F) £ 100, while an upper limit (r) ^ 10 5 
is inferred from the need of having the fireball opaque till 
the acceleration stage is completed. This poorly constrains 
the properties of the outflow, leaving a broad interval for 
physically possible (F). 

In this letter we analyze the kinematic evolution of an 
inhomogeneous fireball, in which the flow is approximated 
by discrete shells with different intrinsic properties (baryon 
load, energy and Lorentz factor) and the outflow energy is 
dissipated through binary shell-shell collisions. We derive 
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general constraints on the relativistic properties of the out- 
flow both in the simplest internal shock scenario, that re- 
quires a single interaction for each couple of shells, and in 
the case of hierarchically developed internal shocks, in which 
successive shell mergers end up in a single (or a few) more 
massive shells, whose interaction with the ISM produces the 
afterglow. The latter scenario is strongly favored by the need 
of dissipating a significant fraction of ordered outflow energy 
(see also Kobayashi et al. 1997). 



2 KINEMATIC OF THE FLOW 

We describe the relativistic outflow as constituted of Ns 
shells, each with its own relativistic Lorentz factor Fi, mass 
Mi and energy E\ = TiMjC 2 . The global properties of the 
flow are described by the averages (F), (M) and (E) and 
their dispersions, by the total duration of the ejection of 
shells T (as measured in the rest frame of the inner engine) 
and by the time interval between the ejection of two succes- 
sive shells At. If the ejection time of a typical shell equals 
the time of quiescent phase between the production of two 
consecutive shells, At — T/2Ns (it is thus assumed that the 
time necessary to eject a fast or a slow shell is the same, see 
however Panaitescu, Spada and Meszaros 1999). 

The average mass, energy and relativistic factor are re- 
lated as 

(M) = jj-rLg ~ 5.6 x 10- 7 E 52 (r)^N^ 2} M , (1) 

where E = 10 52 Es2 erg is the total energy of the outflow Q 
The dynamics of each shell is characterized by an initial 
phase in which it is accelerated to its final Lorentz factor as 
F ~ R/Ro, where Ro is the distance of the shell from the 
center of expansion at the ejection time (see e.g. Piran 1998). 
This phase ends when the final Lorentz factor is reached at 
the 'acceleration' radius: 

Ra = TiR ~ 10 9 r {i , 2} 7? {0i7} cm. (2) 

During the following evolution, the shell coasts with con- 
stant velocity until an interaction with another shell or with 
the ambient medium takes place. Assuming that the acceler- 
ation phase lasts for a time interval negligible with respect to 
the coasting phase, the interaction between two shells with 
Lorentz factors Fi and T2 (> Ti) and initially spaced by a 
time interval At occurs at a radius: 
2 

Ri = 2 a2 T _ 1 T\cAt = 4.0 x 10 13 r{ 12} TiiV { ^ 2} cm (3) 

where the numerical value has been computed for ar = 
r 2 /ri = 2 and assuming At = T/2N S . 

The interaction of the shells with the ambient medium 
is a continuous process and does not happen at a well 
determined radius. However a typical scale can be esti- 
mated as the distance at which the shell Lorentz factor is 
half its initial value, which is reached when the shell has 
swept up an external mass m equal to its rest mass di- 
vided by the initial Lorentz factor (see Meszaros & Rees 

* Here and in the following we parameterize a quantity Q as 
Q = 10 X Q X and adopt CGS units. 



1997). Hence the deceleration radius can be approximated 
as Mi/Ti = (4/3)nR 3 D nm p , i.e.: 

R D ~ 1.2 x lO 16 ^ 3 n~ 1/3 ^~ 2 2yN~ s 1 ^ cm (4) 

where n is the density (assumed uniform) of the ambient 
medium and shells of equal energy have been considered. 
Therefore, the more the shell is relativistic, the smaller the 
deceleration radius is (Eq. while, on the contrary, the 
largest the radius at which internal shocks occur (see Eq. ^J). 
This allows us to put some constraints on the average rela- 
tivistic factor of the flow, since observationally the internal 
shocks set up before the afterglow, i.e. before the develop- 
ment of the external shock. 

To fully describe the kinematic of the flow, a fourth 
transition radius is important, i.e. the radius at which the 
shell becomes transparent (to Thomson scattering). If, in 
fact, the shell becomes optically thin before the acceleration 
phase is completed, the internal energy can escape from the 
shell and the acceleration process is damped. On the other 
hand, if two shells collide when still opaque, they give rise to 
a single shell that is reaccelerated to an intermediate Lorentz 
factor. The transparency radius corresponds to: 



3 COMPARING THE RADII 

Fig. [l] shows a comparison of the four critical radii for 
a typical burst with energy E — 10 52 erg and duration 
T = 10 s. The ISM has been assumed uniform with den- 
sity n — 1 cm -3 . Two cases have been considered: the left 
panel refers to a burst made by a (relatively) small number 
of shells (N s = 100), while for the right panel N s = 25000. 
The smaller value has been derived assuming that each of the 
burst pulses fitted by Norris et al. (1996) is produced by the 
interaction of a shell pair. Norris et al. (1996) find a number 
of pulses N p ^ 50, and hence we have Ns = 2 Np ~ 100. In 
the latter case, instead, the central engine emits shells at reg- 
ular time intervals, equally spaced by the smallest variability 
time-scale observed in GRB light-curves (200/is, Schaefer 
& Walker 1999). In this situation each of the Norris et al. 
(1996) pulses is considered as a blend of sub-pulses, with 
an envelope corresponding to the varying efficiency of the 
central engine: indeed, Walker, Schaefer & Fenimore (1999) 
have shown that millisecond variability on top of the larger 
time-scale modulation is a common feature of GRB light- 
curves. 

In both cases the maximum value of (r) is bound by the 
requirement that internal shocks must happen before exter- 
nal ones, while the minimum value is constrained by the 
transparency condition. This second constraint is less se- 
vere since, as already mentioned, collisions between opaque 
shells would simply cause a reacceleration of the merged 
shell. Note however that this could play a role in prevent- 
ing the presence of very slow shells outside the transparency 
radius R T . It should be also stressed that numerically the 
limit given by the transparency requirement is similar to the 
minimum value F ^ 100, already obtained from the com- 
pactness argument (e.g. Piran et al. 1996). The acceleration 
radius does not impose any significant constraint, showing 
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Figure 1. Value of the four relevant radii as a function of the bulk average Lorentz factor (r) in a burst with energy E = 10 52 erg which 
lasts 10 seconds that produces 100 shells (left panel) and 25000 shells (right panel). The ISM has been assumed uniform with density 
n = 1 cm -3 . The shaded areas correspond to the interval of (r) in which the four radii are correctly nested. If (r) < 35 (100 in the right 
panel) the shells interact before becoming transparent and cannot emit the observed thin spectrum. If instead (T) > 750 (2700 in the 
right panel) the shells are decelerated by the external medium before the mutual interaction occurs, producing an external and not an 
internal shock. 



that from the kinematic and radiative points of view, shells 
could be potentially accelerated even to very high Lorentz 
factors r ~ 10 5 . 

From the condition Ri < Rr> we get an upper limit: 



r Um = 800 7\ ' E, 



3/8 B?V8„-V8 jyl/4 



S,2 



(6) 



This limit is quite robust with respect to E and n while it is 
more dependent of T and Ns ■ A variation of three orders of 
magnitude of the total energy changes it by a factor ~ 3 only. 
The density of the ISM has been assumed uniform and ~ one 
proton cm -3 . This is probably a lower limit. However, even 
in the case of the hypernova scenario (Paczynski 1998), in 
which GRB would occur in a much denser environment (n ~ 
10 4 cm" 3 ), Trim decreases only by a factor of ~ 3. Finally, 
the total burst duration of 10 seconds corresponds to the 
mean T90 parameter of the long GRBs, and is appropriate 
for the majority of bursts. Its small uncertainty influences 
Tjxm rather weakly. 



4 EFFICIENCY OF INTERNAL SHOCKS 

The efficiency of internal shocks in converting the bulk out- 
flow energy into internal energy can be easily estimated since 
each collision satisfies energy and momentum conservation 
(see also Kobayashi et al. 1997). Consider two shells of rest 
masses mi and m 2 and Lorenz factors Ti and r 2 (> Ti), 
respectively. Calling e the internal (random) energy of the 
merged shell after the interaction, we have: 

ri mi + Y2 mi = Ff (mi + m2 + f/c 2 ) 

Ti /3i mi + T 2 Pi m 2 = T / p f (mi + m 2 + e/c) , (7) 

where the subscript / refers to quantities after the inter- 
action and r = (1 — /3 2 )~ 1//2 . If we assume that all of the 



internal energy is converted into radiation, we obtain - in- 
dependently of the emission mechanism - an upper limit for 
the efficiency r\ — e/(Timi + F2m2). From the conserva- 
tion equations an implicit solution for the final bulk Lorentz 
factor of the merged shells and the maximum radiative effi- 
ciency can then be derived: 

Pi + ar a m Pi 



1 + or a m 

r/(l + a m ) 
Fi (1 + a r a m ) ' 



(8) 



where a m = m 2 /mi . The above relations give an upper limit 
to the fraction of energy that can be radiated in photons. In 
fact, strictly speaking, rj is the fraction of bulk kinetic energy 
converted into internal energy. In the standard synchrotron 
shock model, this random energy is equally shared among 
protons, electrons and magnetic field, and only one third of 
this energy (the fraction going to electrons) can be radiated 
(see e.g. Panaitescu et al. 1999) . The major features of Eq. ^ 
are that for a large difference in the Lorentz factors (ar ^> 1) 
the efficiency can approach unity and that for a fixed value 
of ar the maximum efficiency is reached when mi = m 2 . 
For ri,r 2 ^> 1, as for internal shocks, the expression for r\ 
reduces to: 

71 - 1 / ; r ( 9 ) 

y/l + a m (a m + a r + 1/ar) 

An even simpler relation is found if the two shells have the 
same total energy (i.e. ar = l/a m ): 

r] ~ 1 __±+l^ = . (10 ) 
V2 + 2a 2 1 

In this latter case the efficiency is always lower than 30%, 
independently of ar- Fig. ^ shows the efficiency for low- 
intermediate values of ar in two limits: the most efficient 
situation (mi = mi, dashed line) and the equal energy case 
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Figure 2. Radiative efficiency of internal shocks as a function 
of the ratio between the relativistic T factors of the interacting 
shells (Eq. ^ Ti , T2 ^> 1). The solid curve refers to the situation in 
which the energies of the two shells are comparable, i.e. mi Ti = 
m.2 T2 , while the dashed curve shows the highest efficiency for 
a given «r, corresponding to shells of equal mass. Since both 
shells have been assumed highly relativistic, the efficiency does 
not depend on the value of Ti. 




12 3 4 5 6 
Number of interactions 

Figure 3. The simplest scheme for hierarchical internal shock. 
The number of shells required to generate it with N steps is given 
by2 N 



(rimi = V2m2, solid line). The limit 14, r2 ^> 1 has been 
assumed. Note that, in this case, the efficiency does not de- 
pend on the value of Ti. 

Let us now consider the global efficiency of the 
burst. In the 'standard' internal shock scenario (Rees & 
Meszaros 1994) the shells have a distribution of Lorentz fac- 
tors with mean value (r) and width Ar ~ (r) . This in turn 
corresponds to a distribution of ar with (ar) ~ 2 within 
a factor of order unity (which accounts for the particular 
shape of the distribution). This implies (see Fig. |2J) r\ ~ 5%, 
independently of the mass ratio a m . A numerical simula- 
tion for a log-normal distribution of Y with Ar = (F), gives 
r) = 6.8%. 

These (upper limits on) efficiencies are problematic 
since, in the absence of extreme beaming, the total energy 
release of the most powerful GRBs would exceed by some 
orders of magnitude the maximum energy release achievable 
by current models, that involve a stellar mass black hole as 
the energy source. If, e.g., GRB 990123 had an efficiency 
of 5%, the required energy would be £iso — 5 x 10 erg. 
An extremely narrow beaming angle 9 ~ 0.5° would be 
then needed to reduce the energy release to the value Eg = 
10 52 erg. A higher efficiency can be only obtained by allow- 
ing for a broader distribution of relativistic factors F and 
imposing that all shells have roughly the same mass, some- 
how requiring a finely tuned variable engine. Moreover, a 
very broad distribution of Lorentz factors implies that the 
Compton drag effect becomes crucial, as described in the 
following section. 

4.1 Compton drag 

Let us assume that the emission produced in the interaction 
between two shells is isotropic in the shell comoving frame 
K' . Consider also a third shell, faster than the previous two, 
with Lorentz factor F' in that frame. This moves in the radi- 
ation bath of total energy _Ep h produced in the interaction, 



and thus Compton scatters a fraction tt of the photons, in- 
creasing their energy by a factor F' 2 (here tt is the Thomson 
optical depth of the fast shell). The total energy lost by the 
fast shell in the frame K' is hence tt (r' 2 — l)E' pb . 

If tt (r' 2 — 1) approaches unity, there are two major 
consequences: first a significant part of the energy of pho- 
tons is due to the inverse Compton mechanism rather than 
to the internal shock dissipation process; second, the fast 
shell looses a significant fraction of its bulk kinetic energy 
to up-scatter the primary photons and this causes a brak- 
ing of the fast shell (in the frame K'). Thus, even if the 
original distribution of F factors were very broad, the early 
interactions would give rise to an efficient inverse Compton 
drag on the faster shells, reducing the mean value of ar in 
the successive collisions and hence the global efficiency of 
internal shocks. Moreover, in this case the net energy pro- 
duced through inverse Compton scattering would be equal 
or even larger than that directly due to shocks, and the pri- 
mary emission mechanism would be different. We plan to 
investigate this issue in a forthcoming paper. 



5 HIERARCHICAL INTERNAL SHOCKS 

The low efficiency predicted by Eq. ^ refers to a single in- 
teraction between two shells. However, each time this inter- 
action takes place, a new shell is formed which can in turn 
catch up (or be caught by) another one. If the i th interaction 
has efficiency 77;, after N of them the fraction of the initial 
bulk energy converted into photons will be: 

JV JV 

m = ^ ~ *-0-(i - m) - ") <_1 (ii) 

1=1 1=1 

If the efficiency is constant (rji = 77) and TV is arbitrarily 
large, all the energy of the outflow can be converted into 
radiation, as rj^ = 1 for any value of r\. For the typical 
efficiency rj = 6.8% derived above, about 9 interactions are 
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needed to obtain an overall 50% efficiency. Assuming that, 
as sketched in Fig. || the collision 'tree' develops as a binary 
bifurcation, the number N of shells that has to merge into 
a single one after ~9 interactions is J\f = 2 9 = 512. 

This result has two main consequences. Firstly, it is 
unlikely that a burst with high efficiency is produced by 
a small (< 100) number of shells: a scenario in which the 
inner engine emits a large number of shells is thus favored. 
Secondly, in deriving r Un , (Eq. ^ we assumed that a single 
interaction was enough to power a GRB. However, within 
the hierarchical scenario, a shell must interact with N other 
shells and hence Eq. ^| must be modified as: 



: 800 Tf 3/8 EH 8 n~ 1/8 N 1 ^ 4 M~ 1/4 



(12) 



which gives r Um ~ 650 for - say - iV s = 25000 and N = 512. 
Even if this is an illustrative value only (the actual number 
depending on the details of the initial distribution of the 
Lorentz factors), we conclude that values of (F) larger than 
~ 1000 are unlikely. 



6 DISCUSSION 

We have analyzed the kinematic efficiency of internal shocks 
in a fireball made of many individual shells. Since the radius 
at which the internal shocks set up grows with the average 
Lorentz factor of the flow, while the development of an ex- 
ternal shock is favored by a higher relativistic motion, we 
can put limits on the average Lorentz factor (r) of the flow 
if - as inferred from observations - internal shocks develop 
before the external one. Moreover (see also Kobayashi et al. 
1997) we conclude that a simple internal shock scenario in 
which each shell is caught up by a single faster shell, suffers 
from a very low radiative efficiency, unless extremely differ- 
ent Lorentz factors are involved. These would however cause 
a dramatic Compton drag effect, which in turn constrains 
the possibility of a very broad F distribution. A possible so- 
lution is to have a flow with moderately relativistic Lorentz 
factors 100 (r) 600 in which collisions develop until a 
single (or a few) massive shell is formed. For this to hap- 
pen the shells have to be ejected at small time intervals, of 
the order of milliseconds. In this case each burst is made 
of more than a thousand spikes that merge to produce the 
broad peaks often observed in GRB lightcurves. These peaks 
would then reflect a modulation of the hidden central engine 
that powers the burst, rather than a single collision. 

It is interesting to ask whether an outflow with signif- 
icantly higher Lorentz factors produces a burst or not. In 
a flow with very high (F), the internal shocks set-up at a 
larger radius, and hence a hierarchical internal shock cannot 
develop. Since the efficiency is roughly proportional to the 
number of shocks a shell undergoes, a higher bulk gamma is 
linked to a lower efficiency. Observationally, this implies that 
the ratio between burst and afterglow fluences is a function 
of (r) , being larger for smaller (F) . A further consequence is 
that very short (millisecond) bursts should be characterized 
by smaller efficiencies and brighter afterglows. To date, this 
can be only a prediction since the BeppoSAX trigger works 
only for long burst. 

As a final note, following the explosion of the bright 
burst GRB 990123, two values of (r) have been derived. 
Sari & Piran (1999) obtain (r) ~ 200 from the optical flash 



intensity while Liang et al. (1999) estimate (F) ~ 350 from 
the lightcurve properties. Despite the differences, both val- 
ues are in agreement with the general limits derived in this 
letter. 
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